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We use a formal recurrence relation approach to multiple parton scattering to find the complete 
solution to the problem of medium-induced gluon emission from partons propagating in cold nuclear 
matter. The differential bremsstrahlung spectrum, where Landau-Pomeranchuk-Migdal destructive 
interference effects are fully accounted for, is calculated for three different cases: (1) a generalization 
of the incoherent Bertsch-Gunion solution for asymptotic on-shell jets, (2) initial-state energy loss 
of incoming jets that undergo hard scattering and (3) final-state energy loss of jets that emerge out 
of a hard scatter. Our analytic solutions are given as an infinite opacity series, which represents 
a cluster expansion of the sequential multiple scattering. These new solutions allow, for the first 
time, direct comparison between initial- and final-state energy loss in cold nuclei. We demonstrate 
that, contrary to the naive assumption, energy loss in cold nuclear matter can be large. Numerical 
results to first order in opacity show that, in the limit of large jet energies, initial- and final-state 
energy losses exhibit different path length dependences, linear versus quadratic, in contrast to earlier 
findings. In addition, in this asymptotic limit, initial-state energy loss is considerably larger than 
final-state energy loss. These new results have significant implications for heavy ion phenomenology 
in both p+A and A-l-A reactions. 

PACS numbers: 24.85.-|-p; 12.38.Cy; 25.75.-q 



I. INTRODUCTION 

Non-Abelian final-state medium-induced radiative en- 
ergy loss in the quark-gluon plasma (QGP) is the best 
studied many-body perturbative QCD (pQDC) appli- 
cation for high energy nuclear collisions. Several theo- 
retical approaches that address this question have been 
well documented in recent reviews [U, 0, H, 13) Hi- In 
contrast, initial-state energy loss in cold nuclear matter, 
pertinent to hard jet and particle production in proton- 
nucleus (p+A) and nucleus- nucleus (A+A) collisions, has 
not been studied so far. 

Strong motivation for in-depth investigation of cold nu- 
clear matter energy loss comes from the finding that at 
most 1/2 of the forward rapidity suppression in d+Au 
collisions measured at the Relativistic Heavy Ion Collider 
(RHIC) [f^, 13] can be accounted for by leading twist ^] 
or high twist [§] shadowing calculations that take into 
account constraints from deep inelastic scattering (DIS) 
on nuclei. Moreover, at much lower center of mass en- 
ergies at the Super Proton Synchrotron (SPS), a similar 
forward rapidity suppression is established and shown to 
not be compatible with shadowing calculations [To| . 

Energy loss for asymptotic t — ±oo on-shell jets that 
do not undergo hard scattering was discussed in |ll| for 
nuclear matter of extent L ~ tq = 1.2 fm. It has been 
argued that the inclusion of the Landau-Pomeranchuk- 
Migdal (LPM) destructive interference effect [l^, [l^ for 
this regime, in the approximation of infinitely large num- 
ber of soft scatterings, leads to negfigible Ai^rad, 1/3 
of the magnitude of final-state energy loss Final- 
state energy loss in cold nuclei [l| was found to have the 



same qualitative behavior as final-state energy loss in the 
QGP, though much smaller in magnitude and with pos- 
sible relevance to suppressed hadron production in semi- 
inclusive DIS. None of these regimes can yield significant 
and phenomenologically relevant contribution to cold nu- 
clear matter attenuation at collider energies. 

Thus, the primary purpose of this manuscript is to de- 
rive and properly compare the Bertsch-Gunion, initial- 
and final-state energy loss in cold nuclei. The dominant 
contribution to AE^ad can then be used in heavy ion 
QCD phenomenology [31 • The secondary goal of this 
paper is to clarify the prin ciple difference between radia- 
tive and coUisional p^.llB| energy loss in cold/hot nuclear 
matter. Recent studies have made comparisons between 
Ai^rad and AiJcoi, emphasizing a deep LPM cancella- 
tion regime, which is not representative of the process 
of radiative energy loss. In addition, different and often 
incompatible formalisms are used in such comparisons. 

We start by recalling the energy loss results for elec- 
trodynamics (QED) [l3|. Consider medium of atomic 
density p = N/V. Each atom has Z electrons of electric 
charge e and mass m. An incident particle with electric 
charge ze and mass M, E = 7M , p = PjM, undergoes 
multiple Coulomb scattering in such a medium. Its coUi- 
sional energy loss per unit length, including momentum 
transfers down to the mean excitation energy of the elec- 
trons (w), is given by: 



dASc, 
dx 



-^InB, . 



(1) 
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In Eq. m = 27^/3^m/ (cj) and in the high energy limit, 
In Bq ^ (3^, we have neglected a small correction to the 
large logarithm related to the relativistic electron spin. 
The physics behind coUisional energy loss is the energy, 
A£^coi(Q^), transfered to the medium by the incident fast 
particle. Thus, one expects little dependence on the mass 
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of the the incident particle, M, but strong dependence 
on the mass of the target scatterers, m. It is evident 
from Eq. ^ that colhsional energy loss grows at most 
logarithmically with the energy of the incident particle, 
E, and linearly with the size of the medium, L. 

On the other hand, the radiative energy loss per unit 
length is given by 



dl^Erad ^ 1024 3 3 4 . 1 



(2) 



Bertsch-Gunion 






where A — X{M, 7, ujs) may also depend on the screening 
of the target electric charge via the characteristic fre- 
quency ujs- Radiative energy loss arises from the acceler- 
ation of the incident particle, which allows for emission 
of real photons. Thus, one expects that there will be sig- 
nificant dependence on the particle mass, M. Emitted 
photons naturally carry a fraction of the incident parti- 
cle energy, implying that radiative energy loss is propor- 
tional to E or, equivalently, 7. Both features are easily 
seen in Eq. Q and AiJiad grows linearly with L. 

By comparing Eqs. (H]) and ^ one observes that in 
the high energy regime ASrad ^ Ai?coi- Thus, for ultra- 
relativistic particles radiative energy loss is the dominant 
mechanism of momentum attrition. As we will show be- 
low, the same energy dependence of Ai?iad is characteris- 
tic of QCD energy loss. It is only in the deep LPM regime 
for final-state energy loss that the energy dependence of 
AiSj-ad is reduced to logarithmic. Since this is a very 
specific case, care should be taken to evaluate collisional 
energy loss in the same model of momentum transfers 
with the medium to the same power in the expansion in 
1/P+. 

This paper is organized as follows: in Section II we 
review the recursive approach to multiple parton scat- 
tering in both cold and hot nuclear matter, formulated 
in [isl, [1^ . We use final-state energy loss as an exam- 
ple. Next, we derive to all orders in opacity the two 
new solutions for the radiative non-Abelian energy loss 
of incoming partons that may or may not undergo hard 
scattering that produces high-pT or high-i?T particles or 
jets. Section III contains a detailed numerical study to 
first order in opacity of the three different energy loss 
regimes. We identify initial-state energy loss as the dom- 
inant cold nuclear matter AiS^ad contribution, relevant 
to p-fA and A+A collider phenomenology. A summary 
and discussion of our results is given in Section IV. Ap- 
pendixes contain useful kinematic simplifications, which 
allow for the derivation of the all-orders in opacity so- 
lutions given here. Explicit brute force calculations to 
second order in opacity of the new Bertsch-Gunion and 
initial-state Ai^iad are also shown for the purpose of val- 
idating the general results. 



FIG. 1: (Color online) Three distinct cases of medium- 
induced bremsstrahlung are illustrated: 1) Bertsch-Gunion 
case of t — —oo,t = +<x on-shell jets; 2) initial-state energy 
loss in the nucleus followed by a large process, resulting 
in the production of high-pT or high-_ET particles and jets; 3) 
final-state energy loss in the nucleus, after a hard scatter. 



II. RECURSIVE METHOD FOR OPACITY 
EXPANSION 

For radiative processes in QCD, including medium- 
induced bremsstrahlung, it is important to keep track 
of the evolution of the gluon transverse momentum k 
in a plane perpendicular to the direction of jet propaga- 
tion hj. Such a k may arise from single hard or multiple 
soft scattering. The acceleration of the color charges in 
the 2D transverse plane generates color currents whose 
detailed interference pattern determines the strength of 
the non-Abelian Landau-Pomeranchuk-Migdal [l^, 
effect. Let us denote by: 



B,, = H - , (5) 
^(n---»m)(iij2---i„) = ^ ^(jij2---j„) (6) 

the Hard, Cascade, and Bertsch-Gunion propagators in 
the transverse momentum space [H,!!^- In Eqs. ([2]), Q, 
^ and ([6]) k is the transverse momentum of the emit- 
ted gluon and q; are the transverse momentum trans- 
fers from the medium to the jet+gluon system at po- 
sitions Zi. Another important quantity, which enters 
the bremsstrahlung spectrum, is the formation time of 
the gluon, Tf, at the radiation vertex. When com- 
pared to the separation between the scattering centers 
Azj = Zj — Zj-i, which can fluctuate up to the size of the 
medium L, it determines the degree of coherence present 
in the multiple scattering process. We introduce the fol- 
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lowing notation [H, [T^ : 



To = UJQ 



If. 
(ii ■■■I, 



k+ ' 

(k-q.J^ 

fc+ 

(k - q.i - 



k+ 



(7) 
(8) 
(9) 



In Eqs. (O, ^ and ([5]) fc+ is the large positive lightcone 
momentum of the bremsstrahlung gluon. We note that in 
calculating the amplitudes for the gluon emission, prop- 
agators of the type Eq. ([3]) - Eq. © come with a factor 
—2igs£± ■ (• ■ ■) that we don't write explicitly, see Ap- 
pendix El Similarly, at the level of squared amplitudes 

E 



polarizations fJ- 



= — .g^jy- For physical quantities, this 
leaves us with the product of the 2D propagators in the 
amplitude and its conjugate. 



A. Constructing the reaction operator 

In the limit of high energy parton propagation in mat- 
ter, multiple interactions are path ordered. The leading 
nuclear-size-enhanced contribution, oc L/X, to the mod- 
ification of such partonic systems arises from two gluon 
exchanges with the strongly interacting constituents of 
the medium, as shown in Fig. [2l Unitarization of multi- 
ple scattering requires inclusion of three distinct cuts in 
the Feynman diagrams and is also illustrated. We note 
that n > 2 gluon exchanges in the region of local color 
neutralization l/jj, <C A will lead to higher order, oc 
corrections to the in-medium scattering, which are not 
nuclear size-enhanced. These are absorbed in the mean 
free path A for phenomenological applications. It is un- 
derstood that all possible ways of attaching the momen- 
tum exchanges to the constituent partons of the propa- 
gating system should be considered, which increases the 
complexity of the calculation with the number of con- 
stituents. In the simplest case of single parton propaga- 
tion, one can describe its transverse momentum diffusion 
due to the random walk in nuclear matter [2Q|] . Recently, 
the dissociation of heavy mesons in the QGP has been 
calculated by applying this general approach to a qq sys- 
tem [21*1. Finally, the same classes of diagrams, shown in 
Fig. [21 were used to calculate direct photon and dilcpton 
emission to first order in opacity [2^ . 

For the purpose of this paper, we are interested in a 
two-parton, jet-t-gluon, system propagation. The prob- 
lem of radiative energy loss in QCD is more complicated 
than meson dissociation due to the multiple possibili- 
ties for the location of the gluon emission vertex in the 
many-body scattering. We first review [18| the construc- 
tion of the reaction operator, illustrated in Fig. [2l Let 
(x, k; c) be the amplitude when n correlated inter- 
actions between the system and the medium may have 
already occurred. Here, x = /p~^ and k are the kine- 



matic variables and c = Tc is the color matrix associ- 
ated with the radiative gluon. For each scattering center 
the system may: not interact; interact via a single-Born 
scattering; or interact via doublc-Born scattering. These 
correspond to the three possible cuts when an amplitude 
times its conjugate amplitude are considered. Therefore, 
the SU(3) color and kinematic modifications that arise 
from in-medium scattering at position Zm are represented 
by the operators: 

Im — Unit operator (no interaction) , (10) 
Dm — Direct operator (single — Born) , (11) 
Vm — Virtual operator (double — Born) . (12) 

Thus, starting with an initial condition Goix, k; c), which 
can be the vacuum radiation associated with the hard 
scattering of the parton or the lack of such radiation for 
asymptotic on-shell jets [l^, [l^, we can construct the 
amplitude as follows: 



Ai---i„(a;,k, c) = Yl (SoA„A + Sl^i^Drn + 62,i^Vm 



xGoix.k; c) 



(13) 



In Eq. ^i.irr, are the Kronecker symbols and the 

indexes im keep track of which type of interaction, 
Eqs. Uni), pT|) and p^ . has occurred. The operators 
are ordered from right to left as follows: Oi„Oi„_i • • • 
for in > in~i > ••■ > *!■ The conjugate amplitude 
(x, k; c) is then uniquely defined since no interac- 
tion is accompanied by a double-Born term, a single-Born 
interaction is accompanied by a single-Born term, and a 
doublc-Born interaction is accompanied by the unit term, 
see Fig. [H 

Ai-*"(x,k,c) = GS(x,k;c) x 

n 

n {^^,^yl+K^A+52,^J^) ■ (M) 



R = 



l/p, « X 




FIG. 2: (Color online) Representation of the reaction oper- 
ator for an arbitrary propagating system of partons. Three 
different cuts, corresponding to the relevant single-Born and 
double-Born interactions, are shown. 
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In Eq. p4)) the operators are ordered from left to right: 

^li • • • ^L-i^L ^'^^ > > ■■ ■ > h and act to the 
left. 

The contribution to the cross section arising ft'om n 
correlated interactions can then be written as follows: 

Pn{x,k) - ^^-••'"(C)A,....„(C) 
2 2 

= ••• 51 4i...«„(a:^.k;c)Ai...»„(a;,k;c) 

11=0 j„=0 





Col(c) 



in 



a; Col(c) 



qi,ai 



i(COo-Wi) Zj 




(At A, + v; + v;!) Ai...»„-i {x, k; c) . (is) 

From Eq. (jlSp wc can identify the reaction operator, 

i? = (A), i3„ + K + F„t) , (16) 

as the basic input in this iterative approach to multiple 
scattering. 

Clearly, to obtain the contribution from any number, 
n, of correlated interactions to the cross section (or the 
radiative gluon differential distribution) one has to: 

1. Understand the structure of the direct and virtual 
operators, D„ and Vn- The unit operator 1 is triv- 
ial. 

2. Construct the reaction operator, Eq. 

3. Identify the relevant initial condition for the prob- 
lem at hand, oc Glf{x, k; c)Go{x, k; c). 

4. Solve the recurrence relation with this initial con- 
dition. 

We first consider the action of the direct operator l)„ 
at position z„ on an amplitude with n — 1 correlated 
scatterings. In this case, there is a single momentum 
transfer q.^ in the amplitude and the same momentum 
transfer due to (5^(qi— q^) in the conjugate amplitude p^ . 
see Appendix A. The result of such action on Ai-^...i^_i 
can be represented as: 

AiAi---i„_i (a;, k; c) = (a„i + §„ + Ai)Ai---i„_i (x, k; c) 

— Ojn^ii ■■■i-n. — i (.^^ 

- (-1 B„ e^"«^"[c,a„]rei(Ai...^„_J . (17) 

The first term in Eq. (|17p corresponds to a momentum 
exchange with the energetic jet. In the high energy limit 
terms of the type p~ ~ (^ ■ qi)^/p"'" are suppressed and 
we do not keep track of the transverse modification of the 
parent parton. Such interaction corresponds to an addi- 
tional color factor a„. The second term in Eq. ([T7]) arises 
from the momentum transfer q„ to the radiative gluon. 



+ e 

Col(c -> [caj) 

- (-l)\e'«'«^'[cA]Tei 



FIG. 3: Diagrammatic representation of the sum of ampli- 
tudes generated by the direct single-Born scattering DiA. 



which cannot be neglected since the typical k ^ J^i'ii- 
If the gluon emerges with momentum k after the mo- 
mentum transfer, in the amplitude Ai - i„_i it has mo- 
mentum k — q„. The interference phases ^'^(^0-^^^)^^ 
arise from the difference in the longitudinal momentum 
components before and after the the gluon interaction. 
Finally, if the gluon is represented by color matrix c, 
in the amplitude Ai - i„_i color is rotated as follows, 
■j:ca„d^^ [c,a„], i.e. [c, a„]. 

In addition to the above modifications to existing di- 
agrammatic contributions for the amplitude with n — 1 
scattering centers, the acceleration at position z„ acts as 
a source of a new color current, B„, with a phase factor 
e'""^". Thus, the last term in Eq. (|17p represents the 
diagrammatic contributions where the gluon is emitted 
immediately before or after the scattering center. In the 
former case, the gluon may also interact with this center. 
The parent parton has a cumulative color factor 

re/(Ai. J = (an-i)'"-^ • ■ • [aiY' , (18) 
^ („^)2-n . . . (a„_^)2-«„-i . (19) 

In the opacity expansion approach there is a unique re- 
lation between the color factors in the amplitude and its 
conjugate: 



c 



n-1 



(20) 



where Cn is the quadratic Casimir in the fundamental 
or adjoint representations for quark or gluon parent par- 
tons, respectively. We denote by Ny and the number 
of double-Born interactions in the amplitude and its con- 
jugate: 



7V„ =iV„(Ai-j„_J 



n-l 

E 



(21) 
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a; Col(c [c,ai]) 



-J (-Y)\e"^»^-cTe, 



FIG. 4: Diagrammatic representation of the sum of ampli- 
tudes generated by the virtual double-Born scattering ViA. 



n-l 

iV„ = 7V,(^^--'"-0 = 5] <5o,^™ , (22) 

m— 1 

and note that from a multinomial decomposition of zero 
we have: 

ii , ■ ■■■,im 

= (-i-i + ir = 0. (23) 

The factor (—1/2)^" arises from the symmetry in the two 
gluon legs at the same location Zm and from the fact that 
when both momentum exchanges are in the amplitude or 
its conjugate we have = (— i)^ = —1. 

To summarize, for the medium-induced radiative prob- 
lem in QCD, the single-Born or direct interaction can be 
represented as follows: 

Dn^{aJ + Sn+Rn) . (24) 

Here, 

S„ = X e'^•^o-^^)z„^^q„■b ^ ^25) 

with Id = i Vi^ and e'i-^/(k) = /(k - q). It is im- 
plicit that the color rotation if'^°'"''-d, yielding [c, a„], is 
of the appropriate color matrix d not shown explicitly in 
Eq. ([^5]) . The additional Bertsch-Gunion operator reads: 

xB„e'-''^" Tei(Ai...z„_J . (26) 

Next, we consider the double-Born interaction of the 
jet at position z„, as shown in Fig. [H In this case 



there are two equal and opposite momentum transfers 
Qi in the amplitude (or the conjugate amplitude) due to 
(^^(qi +cti) [la], see Appendix A. The modification of the 
amplitude (a;, k; c) is found to be: 

V^iA^j^ . . . (x, k, c) — ^ yl^^ . . (x, k, c) 

_ (_i )^"(-^----i) ^B„e-«^"cTe,(A,...,„_J . (27) 

The first term in Eq. ([?f|) corresponds to the case when 
both momentum exchanges are with the parent parton or 
the radiative gluon. In this case there is no net momen- 
tum transfer, no additional phase factors arise and the 
two color matrices yield the quadratic Casimirs Cr and 
Ca- The factor (—1/2) was discussed above. The sec- 
ond term in Eq. (|27p represents the situation where one 
of the gluon legs is attached to the jet and one to the 
bremsstrahlung gluon. Since there are two possible com- 
binations 2 X (—1/2) = —1 . The shift in the transverse 
momentum k k — q„, the phase factor e*'"""""'^" 
and the color rotation if^°-^'^d, i.e. c — » [c, a„], are all 
the same as in the second term in Eq. (fT7|) . The differ- 
ence is that the interaction with the parent parton gives 
an additional color factor a„. Finally, parts of the di- 
agrams where a gluon is emitted immediately before or 
after the interaction point z„ are combined in the last 
term in Eq. ((27|) . This term is identical to the last term 
in Eq. (fT7|) . except for the color factor since we have car- 
ried out the simplification a„[c, a^,] = — (Cyi/2)c. The 
additional (—1) arises from the last two gluon exchanges 
at z„. In summary, the double-Born interaction at Zn 
can be implemented by the following operator 

Vn = —^{Ca + Cjl)l ~ anSn — O-nBri 

= -ar,Dn - HCa - Cr) , (28) 

and is clearly not linearly independent of the single-Born 
term, Eq. (p4)) . 

Now, we can proceed and construct the reaction op- 
erator Rn to relate the n"^ order in opacity gluon emis- 
sion "probability" distribution Pn{x, k) to the probability 
at (n — I)*'* order. We express the reaction operator in 
Eq. (fT6|) as follows: 

Rn = (Dn - an)HDn - an) - CaI 
= (Sn+Bn)HSn+Bn)-CAi 

= (SiSn -CAi) + BiBn + 27^e SlBn . (29) 

Noting that i/ca„d(_j) jca„d' ^ CASd,d', we see that the 
first term in Eq. 

{S^Sn - CaI) = CA{e-''-^-e-'^-^- - 1) , (30) 
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gives a homogeneous contribution to the functional re- 
currence relation: 

^^-■■^""1 {x, k; c){SiSn - Ca)A,...^^_, {x, k; c) 
= Ca {Pn~i{x,k- q„) - P„_i(.T,k)) 

= C^(e'i-^-l)p„_i(a;,k) . (31) 

The second term does not contribute beyond first order 
in opacity, n ~ I, since 



^^-■■^"-1 (x, k; c)BtB„Ai. (x, k; c) = |B„|2 



(32) 



For this proof we used Eq. ((20|) . and the identity ([23l) . 
For n — 1, however, the CijC^i |Bip contribution to 
^ri-i(2;,k) survives. Finally, the non-diagonal term in 
Eq. dim reads: 

27^e ^i-"-! {x, k; c)Bt 5„Ai...z„_i (i^, k; c) 

= -2Ca B„ • (Re e-^""^"e*'i"-^I„_i) , (33) 

where we have used Eqs. and Writing explicitly 
the expression for I„„i and representing the (n — I)*'' 
step in this iteration via the Sn-i and Bn-i operators 
we obtain: 



- E (-i)""^-'*-""-^^J^(A....._) 



,i(wo-w„_l)z„_i iq„_i-b 



l)c+ [a„_i,c]B„^i 



In Eq. p4p the term proportional to Bn-iAi-^...i^_2 van- 
ishes for n > 2. The proof, up to an insignificant differ- 
ence in color factors, is again based on the multinomial 
decomposition of zero, Eq. ([^5]) . 

From Eqs. ([50)1 and ([55]) . the basic iteration step for 
the full solution of the problem of medium induced gluon 
radiation can be written as follows: 

P„(x,k) = CA(P„-i(x,k-q„) -P„-i(2:,k)) 

-2CaB„ • (jZe e-^""^"e''i" *'I„_i) 

+Sn,iCACn]Bi\^ ■ (35) 
The inhomogeneous term in Eq. ((35)) is expressed as: 



In-l — 



CA^e 



!(wo- 



-1> 



„iq„_ib 



-I I 



'-n-2 



(36) 



(34) 



^ -(5„_iaCACflBie^'^°^i , 
which is a direct consequence of Eq. p4p . 

We emphasize that each radiative problem yields a 
different solution, related to the boundary condition 
G'o(x,k;c). For final-state radiation this initial con- 
dition is given by the bremsstrahlung associated with 
the hard scattering of incoming partons. This case has 
been considered in detail in Ref. [l^ as the first com- 
plete application of the GLV approach. Explicit solu- 
tion to all orders in the correlation between the momen- 
tum transfers from the multiple scattering centers to the 
jet-|-gluon system, suitable for further analytic approxi- 
mations or numerical simulations, has been found for the 
case of final-state radiation, the third case illustrated in 
Fig. [TJ We emphasize that these are not correlations be- 
tween the scattering centers of the target arising from the 
strong nuclear force. Subject to the coherence criterion, 
Tf > max(|zi — Zj\,0 < i,j < n) at n*** order in opac- 
ity, the momentum transfers are correlated in the sense 
that the bremsstrahlung from the hard collision (at po- 
sition zq) and the multiple soft interactions (at positions 
Zj > zq) interfere to give a contribution to the medium- 
induced spectrum of gluons. For completeness, we quote 
the result here llSl: 



_ dN9{FS) _ Cro, 
dk+d^k ~ 7r2 



■E 

n=l 

2C(i 



dAz 



1 dael{Zj) ^2 



(!,...,„) • 2^ B 

(m+l,---,n)(m,---,n) I ^OS 
m=l \ \ k=2 



l{zj) d'^qj 

^^{k,-,n)^Zk -cos E 

Azk 



,(37) 



where ^2 = a-i^d B(„^x „) = B„ is understood. In the case of final-state interactions, zo ~ is the point of the initial 
hard scatter and zl — L is the extent of the medium. The path ordering of the interaction points, zl > Zj+i > zj > zq, 
leads to the constraint X]"=i '-^^i — ^l- One implementation of this condition would be Azi G [0,zl — ] 
and it is implicit in Eq. p7p . If we write the 2D propagators and interference phases in Eq. p7p more explicitly, we 
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have: 



dN3{FS) _ Cna,^ 
k ^^^^ - —^2^ 



dk+dP'k 



n=l 



n p n ^ / 

J 7 = 1 



k q 



It 

1 — • • ■ — 



k q 



m+l 



k q„ 



(k-qi q„)^ ^iV(k-qrn+i Qn)^ (k-q,„ q„)^ 



X COS 



E 

\k=2 



(k - q/c 



-^Azj j - cos 1^^^ 



(k - q/c 



fc+ 



(38) 



Jet quenching in the QGP, the best known apphcation of The solution for the In-i part in the inhomogeneous term 
non-Abehan energy loss in heavy ion reactions, is based of Eq. (|35p reads: 
on such a solution. For details, see 



B. Bertsch-Gunion radiation to all orders in 
opacity 

The first new solution that we obtain is this paper is for 
the case of asymptotic on-shell jets, initially considered 
by Bertsch and Gunion flTl and illustrated as the first 
case in Fig. [TJ Although it is not directly applicable to 
the physics situation of high-pT particle production due 
to the lack of hard scattering, it is a necessary step to 
fully solve the problem of initial-state energy loss. The 
absence of hard bremsstrahlung yields a simple initial 
condition: 

Go(a;,k;c) = . (39) 



10 = 0, (40) 

11 = -CACflBie'"""^ , (41) 

I„_i = -CrC^-^ Yl (e*('^°-"'")^"'e*i'"-^ - l) 

m=2 

X Bie''"°^i . (42) 



Consequently, the solution for Pnix, k) can be written in 
the form 



Pnix,k) = CrC:1 n (^''''"■^ - l)|Bip + 27^eCflC^ ^ 



m=2 



m=2 



Jl (e'l^ b-i) 



m — 1 



(43) 



J=2 



In Eq. (f43|) we recall that, for the special case of ti = 1, Yl^i=2 — ^ ^'^d 01=2 — ^- ^^^y verify that our result 
is a solution of the master recurrence relation Eq. (P5|) by rewriting it in the form: 



P„(x,k) = CA{e'^--^-l) 



ri-l 



m—1 



X \Y C'A(e*^"°~'^^)^^"e*'i^-'" - l)Bie''^«^i 
X Bie*"°^i . 



n-1 



B„e*i'"-''e-*'^''^" 



i=2 



(44) 



r 



To obtain the contribution of n correlated scatterings age over the momentum transfers q„ in Eq. (|43p . Let 
to medium-induced gluon production we have to aver- 
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{\/aei)daei/<P'<\m be the differential distribution of mo- 
mentum transfers at position z^. We note that the first 
term in Eq. (|43|) can be written as follows, 



i=l 



n / 



1 da, 



]-[(e»q™-b_i)|B^| 



m=2 



.1—1 

n 



(45) 



m=2 



In Eq. (j45p we are able to carry out the simplification 
including the n = 1 term since 



(46) 



The inhomogeneous term in Eq. (|43|) can also be sim- 
plified. For the momentum transfers qm+i, • ■ ■ , Qn the 



result follows from Eq. p5)) . For the momentum trans- 
fers q2, • ■ • , qm-i we use 



m— 1 



n / ^'q. 

"m — 1 



1 d(T, 



m — 1 

i=2 



I da el ^2 1 \ 

(5 (q») 



1=2 
m— 1 



i=2 



(47) 



We note that the terms n — 1 and n ^ m can also be 
included in the general representation, Eq. ([T7)) . due to 
Eq. (|46p . Thus, we write the contribution of the multiple 
scattering centers as follows: 



J 



n / ^^q^ 



1 dap 



a el d^q, 
-t-27ee ^ 



n / '^'q^ 



1 daei 
CTel d^qi 



Y[ e"!"^ Bi 



m=2 



n ' 

j—rn+l 



i=2 



(48) 



We note again that in Eq. (|48p the products of momen- 
tum shift operators or phases and momentum shift op- 
erators are applied sequentially from the left with each 
operator acting on the function resulting from the previ- 
ous one. 

The last step in writing the explicit solution is to 
carry out the action of the momentum shift operators 
in Eq. The first homogeneous term can be easily 

simplified since 



X e ^j=2 
= B 



B 



(2---m-l)(l---m-l) 



(m+l---n)(m---n) ' ■'-»(2---n)(l---n) 



xe 



(50) 



W e"!"'*' |Bip |B(2. ..„)(!...„) 



(49) 



The inhomogeneous term can also be simplified as fol- 
lows: 

n 

g«q.-bB^g«q™-bg-«u;oz„ 

j—va+l 

m — 1 

i=2 



B 



{m+l---n)(m---n) 



e 



The overall normalization for the differential gluon dis- 
tribution is set by the color factor, the strong coupling 
constant and a phase space factor that combine to pro- 
duce the factor CflOfs/Tr^. We note that Cr signifies the 
color charge dependence of the gluon bremsstrahlung in 
the small energy loss limit AE,-g^^/E <C 1. Similarly to 
the case of final state gluon bremsstrahlung, in the limit 
of small lightcone momentum fractions x and small trans- 
verse momenta k the result is "color trivial" , retaining 
only the quadratic Casimirs in the adjoint representa- 
tion. These enter the mean free path of the propagating 
jet-|-gluon system, indicating that only gluon rescattering 
is important, i.e. Xg{zi). Putting everything together we 
find: 
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dk+d^k 



n 



dAz 



1 



ei{zj) d'^qj 



X B(2. ..„)(!...„) 



B(2...„)(i...„) + 2^ B(j+i...„)(j...„) COS y^cj(j-...„)A2 



(51) 



A direct comparison of this general result to the brute force calculation up to second order in opacity can be found in 
Appendix IB] The integration limits on the separation between the multiple scattering centers were discussed in the 
previous section. Since in this case there is no hard scattering, the maximum separation corresponds to the physical 
size of the medium. More explicitly, our result reads: 



_ dNs{BG) ^ CrOs 
dk+d^k 7r2 



Y[ / dAz, (7^i{zi)p{z^) Jl / d^q_j i- 

,i=l J .7 = 1 ^ ' 



1 da,,(z,)_^2(q^.) 



k - q2 



k qi 



(k-q2 



q„)2 (k - qi 



qn)^ 



i{zj) d'^qj 
k-q2 



k - qi 



-2E 

1=2 



k qM 

(k-q 



i+l 



k - q, 
(k - q^ - 



qn 



(k - q2 q„ 

^ (k_q, 

,J=2 



|2 (k - qi q„)2 

qn)' 



-Azi 



(52) 



We recall the convention B 



(n+l,n) 



B„. Thus, to first 



order in opacity the Bertsch-Gunion case of asymptotic 
t = — oo to t = +00 on-shell jets yields a gluon radiative 
spectrum oc |Bip and no coherence effects. Note that in 
Eq. ((52|) (Tei is the gluon scattering cross section. 



C. Initial state radiation to all orders in opacity 

For hadronic reactions where high-p^/high-ii^T parti- 
cles/jets are detected, the relevant initial-state interac- 
tions that lead to energy loss are illustrated as the second 
case in Fig. ([T]). The difference from the Bertsch-Gunion 
case is that there is always radiation associated with the 
hard scatter at position zl- In particular, it can be writ- 
ten as a boundary contribution in the absence of soft 
momentum transfers from the medium in the form: 



H{x,k;c) = +He^' 



(53) 



Such a term will always be present in the opacity expan- 
sion of the amplitude and its conjugate but modulated 
by the color and symmetry factors associated with the 
preceding interactions of the parent jet with the nuclear 
matter. More specifically, 

Ai---j„(a;,k;c)/5 = (a;, k; c)sg 

■ -, N iv„(^.i....„) 

crez(A,....JHe'"°^^ , (54) 



(^x , k, Cj 



IS 



A' 



'(a;,k; c) 



Tj;(Ai-'")He- 



(55) 



The contribution to the differential distribution at n^^ 
order in opacity is proportional to 

Pn{x,k)is = A'^"''"{x,k;c)isAi^...i„{x,k;c)is 

1 2 

= Tr ^ • • • ^ ^J^...,^(a;,k;c)7sAii...j„(a;,k;c)/s 

ii=0 j„=0 



^i'"*" {x, k; c)_BGAi---j„ {x, k; c)bg 

2 2 . 

H^TrE--ER 



ii=0 i„=0 



2 Tie 



xTj,(A--'")ccrei(-4.l-.J 

2 2 / -,^^V„(^*: 

2, 



He-~Tr^...^ [-\ 



ii=0 i„=0 



xT],(^^i-'")c A,....„(x,k;c)BG 



(56) 



The first term in Eq. (f56|) is exactly the same as for 
the case of asymptotic jets. For this, we directly use 
our results from the previous subsection. The second 
term in Eq. (|56p cancels exactly at any order in opacity 
n > 1. The proof is again related to the multinomial 
decomposition of zero: 



2 2 . 

H^TrE--ER 



BG 



ii=0 i„=0 
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xTll{A''■■■'-)ccT,l{A^,...^J 



(61) 



= C^ixQH2(-^-^ + l 



1 1 

2^2 



. 



(57) 



The last term in Eq. ([56|) is the one that differentiates 
the case of initial-state medium-induced radiation from 
the case of fully on-shell jets. We represent this term as 
follows: 



2H- 7^e e" 



J„ ) 



(58) 



where, similar to the case of I„_i from Eq. ([M]) . the 
expression for J„ can be simplified as follows: 



xc^ii...i„(x,k; c)bg 

n ■■■in-i 
X ^ii...i„_i(x,k;c)sG • 



(59) 



As in the case of Bertsch-Gunion radiation, the terms 
arising from _B„ cancel for n > 1, see Eqs. ([32]) and (p4|) . 
The initial conditions for the iterative solution are the 
Bertsch-Gunion terms that come from the direct and vir- 
tual contributions: 

Ji = -Bie''^°^i aic[c,a{\ - ^Bie^'^"^! cc 

= ^CRCA^ie'^"'" ■ (60) 
The solution for J„ can then be expressed as follows: 



i=2 



(uJo- 
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Substituting Eq. (pT]) in Eq. ([55)1 and recalling the 
integration over the distribution of the momentum 
transfers from the medium, we simplify the new non- 
homogeneous term in the solution for the medium- 
induced bremsstrahlung: 



n „ 

n / d"^^ 



1 da, 



2 CrCI H • TZee- 



i=2 



1 dcTel 



J=2 



1 rfcTei 



- 1 



el rf^q 

xH•B(2...„)(l...„)7^ee"'^^=2 '^.-^A^, 



(62) 



Here, we have denoted Zn+i = zl and [.j„+i_„ = ljo- 
With this solution for the inhomogeneous term, we can 
now write the solution for initial state energy loss as the 
sum of the Bertsch-Gunion case and the destructive in- 
terference term, note the "— " sign in Eq. (p^ . The net 
result reads: 



k+ 



dN3{IS) _ Cnas 



n 



dAz 



\J Xg{Zr) 



JJ / d\ 



1 da^^iz,) 



aei{zj) d^c\j 



|B(2...„)(i...„)P -I- 2 B(2. ..„)(!...„) • ^B(j+i...„)(i...„) COS y^a;(j...„)Azj 



-2H • B(2. ..„)(!...„) COS I ^ w(J■...„)AzJ 
.J=2 



(63) 



To second order in opacity, direct comparison is made in Appendix [Cl with the brute force calculation. For the case 
of initial-state energy loss, zg is the position where the jet enters the medium and zl is the position of the hard 
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interaction. Writing the color current propagators and interference phases directly, we obtain: 



_dN9{IS) _ Cnas 
dk+(Fk ^ 7r2 



.1=1-^ J 7 = 1 



k-q2 



■E 

i=2 



(k - q2 - • • 



(k - q, 



+1 



k - q2 



k2 V(k-q2 



k - qi 



q„)2 (k - qi 



k qj 



q„)2 (k - qi - 
q„ k qi 



(k-qi 



1 daeK^.)_^2(q^.) 



q,i 



qn)2 

• ■ q 



k - q2 



k qi 



(k-q2 



q„)^ (k - qi q„ 



qn)2 

• - qr: 



qn)^ 



E 

ri+1 

E 



(k - qj 



qnj 



fc+ 



(k - ^J 



k+ 



-Azi 



r 



(64) 



III. NUMERICAL RESULTS 

In this section we carry out numerical simulations of 
the different radiative energy loss regimes, Eqs. ([57)) . ([5T|) 
and (j63[) , to first order in opacity. The momentum trans- 
fers from the medium are given by the 2 — > 2 differential 
scattering cross section calculated in the Born approxi- 
mation using a finite range, ~ ^/^J■^ Yukawa potential: 



d(Jei _ CrCt 4a;2 



1 



dA TT (q2+^2)2 



(65) 



Here, Cr and Ct are the quadratic Casimirs of the jet 
and target representations and dA is the dimension of the 
adjoint representation. For QCD, 



CrCt _ 9 1 2 
dA 8 ' 2 ' 9 



(66) 



for gg gg, q{q)g q{q)g and q{q)q{q) qiq)qiq), 
respectively. 

The total scattering cross section a^i can be absorbed 
in the mean free path, Ag = 1/aeiP, and the normal- 



ized momentum transfer distribution from the medium 
is given by: 



1 da, 



aei d^q 7r(q2 -|- ^2)2 



2 

fr X — 



)2 

J max 



)2 

- max 



(67) 

67|l , /i2jj arises from the finite range of integration 



In Eq 

q^ < Qmax — s/A — fiEjct/^. This constraint is relevant 
for small jet energies. 

Next, use the results to first order in opacity derived 
here. The reason for this n — 1 choice is twofold: firstly, 
the formulas are simple enough to allow analytic in- 
sight in the different behavior of initial- and final-state 
A£^i-ad- Secondly, the n = 1 Bertsch-Gunion result is 
the prototypical medium-induced radiative energy loss 
in QCD which sets the scale relative to which the de- 
structive LPM interference effects have to be evaluated. 
For final-state energy loss it has also been demonstrated 
numerically that the the sum of all contributions up to 
3'''^ order in opacity is well approximated by the n — 1 
term [T9|. From Eqs. ([51]), §3^) and jH 

in a medium of constant density we have: 



dk+d'^k 7r2 Jg Ag Jq^<s/A crei d'^q ^ 

\J^^<s/i ^7r(q2+^2)2 [k2(k-q)2j ' 

,.dN3(IS) CrUs dAz f ,o 1 ddel nx. ,9 X, .... ^^^ 

fc+ ,,+ ,J = / ^ / Bip-2H-Bicos c.o i-Az 

dk+d^k 7r2 Ag Jq^<s/i (Jei a^q 

^ CrO^ f 2 \L q2 ^ q2 - k ■ q fc+ 

Jc^<s/4 ^4q^+/i^)HA,k2(k-q)2 k2(k-q)2k2Ag 

^ ^ — — \ \ ^ h2Ci -Bi 1 -cos wiAz 



12 



q2<s/4 



7r(q2 + ^2)2 



2k q 

k2(k-q)2 



L 



(k 



_fc+ 



■ sm 



(k-q)^L 

fc+ 



(70) 



Let us examine qualitatively the behavior of the radia- 
tive spectrum and the total lightcone momentum loss, 
Ap+. For energetic jets {Ejct/mq oo) we use "energy 
loss" and "lightcone momentum loss" interchangeably. 
The magnitude of Ap+ depends on the available phase 
space for the radiative gluon at large . For the case of 
Bcrtsch-Gunion radiation, Eq. k'^ does not appear 

in the integrand. Therefore, in addition to the linear de- 
pendence on the path length L, the e nerg y loss is propor- 
tional to the available energy. In Ref. [ll'l this can be seen 
through the flat rapidity dependence of dN^ (BG) / dyd^k, 
since y = lnA;+/p+. The result is qualitatively similar to 
bremsstrahlung energy loss in QED, see Eq. ([2]), and we 
find: 

Ap+{BG) ~ CRa,p+^ X /bg(a*, Qo,P+) • (71) 

In Eq. (|7ip /bg(MjQoiP^) depends on the momen- 
tum transfers from the medium, the mass scale in the 
monopole nucleon form- factor, Qq, and may have addi- 
tional weak logarithmic dependence on p+ via Q^^^. 

The qualitative behavior of initial-state energy loss can 
be understood in the limit of large p+ when k^i/fc+ <^ 1. 
Expanding the sin(- • ■) term in Eq. (j69p . we find 

Ap+ilS) ^ CRasP+^ X //s(Ai,Qo,P+) , (72) 

with fis{tJ;QQ,P^) < /sg(m, <3o,P+)- The difference in 
the overall magnitude of Ap'*" arises from the cancella- 
tion of part of the Bertsch-Gunion radiation through its 
interference with the bremsstrahlung from the hard colli- 
sion. Still, in the asymptotic regime AE'iad remains pro- 
portional to the jet energy and depends linearly on the 
size of the nuclear matter. This result is qualitatively 
and quantitatively different from the argument given in 
Ref. [3] for the Bertsch-Gunion case in the limit of a very 
large number of soft scatterings, namely that initial-state 
Ai?rad is independent of the energy, grows quadratically 
with L and is smaller than final-state energy loss by a 
factor of three. The reason for this difference is twofold: 
firstly, as indicated above initial-state energy loss was 
approximated by the case of on-shell jets ±oo in @. Sec- 
ondly, while our approach is general enough to address 
both L/\g ~ few and L/Xg 3> few cases, the numerical 
results presented here address the limit of few scatterings 
which we anticipate is relevant for finite nuclei. Whether 
the distinction between initial- and final-state energy loss 
will be reduced for very large number of jet-medium in- 
teractions is still an open question. 

Final-state non-Abelian energy loss has been analyzed 
in detail in (Tsl . [l^ . Naively, one may expect that in the 
high energy regime the expansion of the sin(- ■ •) term in 



the integrand of Eq. ([70]) may lead to an ~ depen- 
dence of the radiative AE'rad- However, when properly 
weighted by the available phase space for the emitted 
gluon, this is reduced to a quadratic dependence, at most. 
In addition, the extra powers of l/k'^ cancel the linear 
dependence on the of the energy loss on i?jet, leaving only 
a logarithmic dependence. Thus, 

Ap+(F5) ~ Cnas^ x /i.s(M, Qo,P+) ■ (73) 

To study the radiative energy loss quantitatively, we 
have to identify the |k| scale, at which non-perturbative 
effects can become important. For cold nuclear matter, 
in the original work of Bertsch and Gunion [ll[, this 
scale was approximated by Qq — nip ~ 770 MeV. In 
generalized vector dominance models, the effective mass 
is bigger, Qo — 1 GeV, due to the contribution of the 
heavier vector meson states. One can also understand 
the appearance of form-factors in a purely partonic pic- 
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FIG. 5: (Color online) Fractional energy loss for massless 
quark jets versus the jet energy, -Ejot, in cold nuclear matter 
of length L — 5 fm. Two different sets of typical momen- 
tum transfer per scattering and gluon mean free path have 
been used for comparison: {fi,\g) — (0.35 GeV, 1 fm) and 
ill, Xg) = (0.7 GeV, 4 fm). 



ture. As |k| 0, the transverse size of the radiative 
gluon A± = 1/k^ exceeds the target parton size and the 
strength of the interaction of the jet + gluon system is 
hmited. For cold nuclear matter 1/Qo 
We take Qq — rriN — 0.94 GeV from the inde pen dent cal- 
culation of dynamical nuclear shadowing [9|,Iio[, sugges- 



Bertsch-Gunion E-loss 



Initial slate E-loss 



Final state E-loss 



< 1/^1 < 1/Aqcd- 
ade pen d' 

five of partonic spots of size ~ 0.2 fm inside the nucleon. 
We implement this mass scale as follows: 



(74) 



The last term in Eq. ()74|) arises for heavy quarks and 
was discussed in [23l |. In addition, at any fixed order the 
opacity expansion series is not a perfect square of a sum 
of amplitudes. Still, the integrands in Eqs. ((69|) and (iTOl) 
can be represented as (|Bip — correction). The require- 
ment that the interference does not cancel differentially 
more than the available bremsstrahlung, induced by the 
medium in the absence of coherence, can be formulated 
as follows to any order in opacity: 



E 

1=1 



k+ 



dk+d^k 



> 



(75) 



Eq. ([75]) can also used to identify and eliminate the parts 
of phase space where the approximations that we made 
are least reliable (e.g., the part of phase space where a 
gluon emitted in the direction opposite to the momentum 
transfer from the medium). To set the kinematic limits, 
we require that the positive gluon energy > Qq and 
the gluon rapidity, 



1 k+ 1 (fc+)2 



be within the rapidity gap from the target to the pro- 
jectile, < Ug < In (pp^.^j^^^^/TOAr). Such kinematic con- 
straints correspond to 



Qo , k+ 

— r < a; — - < 1 

xrriN < Ikl < fc^ 



(77) 
(78) 



We have checked that in the limit of large p"*" there is 
little or no sensitivity of the energy loss to a factor of 
two increase or decrease of the lower bounds in Eqs. ((77|) 
and (i75|) . 

In Fig. [5] we show numerical results for the fractional 
lightcone momentum loss Ap+/p+ of massless quark jets 
in cold nuclear matter of length L = 5 fm. For the refer- 
ence Bertsch-Gunion case, Eq. (|55)) . in the limit of large 
jet energies we recover Ap+ cx . We study two sets 
of momentum transfers and mean free paths describing 
cold nuclear matter: (/ii — 0.35 GeV, A^^ i = 1 fm) and 
(M2 — 0.7 GeV, Xg^2 —4: fm). Glearly, the energy losses 
for these cases differ by close to a factor of two, indicating 
that (/i^/Ag)i = {n^ /\g)2 is not a universal parameter re- 
sponsible for the quenching of jets. It should be noted 
that for only a few scatterings in the medium with small 




FIG. 6: (Color online) Fractional energy loss for massless 
quark jets of Bjct = 20 GeV, 100 GeV and 2000 GeV. Three 
different path lengths L = 2 fm, 5 fm and 10 fm were 
used. Power law fits indicate the system size dependence of 
Ap+/p+ = AE,^d/E. 



momentum transfer (x fi — 0.12 GeV^ quarks can lose 
on the order of 30% of their energy. 

The physically relevant case of initial-state energy loss, 
Eq. is shown in the middle panel of Fig.O Destruc- 
tive interference effects can lead to a large six-fold reduc- 
(76) tion of the fractional energy loss Ap+/p+ for p+ — > oo. 



Still, qualitatively the behavior is similar to that of the 
Bertsch-Gunion case discussed above. For medium pa- 
rameters {fii = 0.35 GeV^ , Ag, 1 = 1 fm), directly 
relevant to phenomenology [ij], quark jets lose ^ 5% 
of their energy. In contrast, for final-state energy loss, 
shown in the bottom panel of Fig. (U the cancellation of 
the radiation in the asymptotic domain is still effective, 
leading to Aj)+/p+ 0. Although there is a 30% dif- 
ference between the two calculations, only in this regime 
of final-state radiation in cold nuclear matter may fJ.^ /Xg 
be considered an approximately relevant jet quenching 
parameter. 

Numerical results for the path length dependence of 
jet energy loss in the three different regimes, Eqs. ((68)) . 

and ([70]) . are presented in Fig. [S] for three jet en- 
ergies, Ejct = 20, 100, 2000 GeV, respectively. By def- 
inition, the reference Bertsch-Gunion case without in- 
terference yields a linear dependence of Ap+/p+ on the 
path length L. Numerical calculations were carried out 
for three different path lengths, L = 2, 5, 10 fm, and 
power law fits were used to extract the power index n of 
the Ap+ oc I/" dependence. At large p'^ initial-state en- 
ergy loss approaches a linear behavior, n = 1, while final 
state energy loss shows approximately quadratic behav- 
ior, n = 2, for static nuclear matter. 



L4 



Both results have important implications for heavy ion 
phenomenology. Let us define the nuclear modification 
factor, 



Rab{pt) 



Tab da^p/dyd^pT 



(79) 



which is used to identify dense matter effects on par- 
ton propagation. In Eq. (jTO]) . Tab is the nuclear over- 
lap function, which relates the hadron multiplicity in the 
many-body collision to the corresponding cross section in 
p-|-p reactions. In A-l-A collisions, the quadratic depen- 
dence of final-state energy loss is reduced to linear, when 
Bjorken expansion is taken into account, and we predict 
the following suppression pattern: 



In Raa{pt) 



Ar2/3 
-KAA^^pa.xi 



(80) 



In Eq. ([50]) iVpart is the total number of participants and 
hab is a microscopic coefficient, which depends on the 
properties of the dense latter. Initial-state energy loss, 
on the other hand, approaches linear dependence on the 
system size TVpait (target), which can be used to predict 
analytically (14ii] the centrality dependence of forward ra- 
pidity hadron suppression in p+A reactions [10]. Here 
A^part (target) is the (average) number of nuclear target 
participants along the projectile line for a given (central- 
ity class) impact parameter. The quantitatively different 
behavior of initial-state energy loss leads to the following 
result in p-l-A reactions: 



\tlRpa{pt) = -KpA Apart (target) 



(81) 



Before we proceed to more differential bremsstrahlung 
distributions it is important to emphasize that the 
quadratic dependence of final-state energy loss on the 
size of the static medium, oc L^, does not imply that 
the magnitude of the energy loss itself is large. On the 
contrary, such dependence arises from the maximally ef- 
ficient cancellation of the large x = /p'^ radiation and 
for the same parameters of the medium, in the limit of 
p^ — > oo, final-state energy loss is negligible when com- 
pared to initial-state energy loss. 

One can gain insight in the energy dependence of 
Ap+ by investigating the radiative intensity spectrum 
xdN^/dx, obtained by integrating the double differen- 
tial distributions Eqs. ^E^, and (HH) over k and 
shown in Fig [T] We have chosen i = 5 fm (/ii = 
0.35 GeV, Xg^i — 1 fm) and the same jet energies 
as in Fig. \6\ For the Bcrtsch-Gunion case there is an 
approximately flat, up to the edge of phase space con- 
straints, dependence of the bremsstrahlung intensity on 
X, leading to Ap+/p+ « constsc. For initial-state energy 
loss there is a finite cancellation of the intensity spec- 
trum at large values of fc+, which leads to a suppressed 
Ap+Zp"*" « const/5. For final-state energy loss, the spec- 
trum is progressively more suppressed by destructive in- 
terference effects that lead to a power behavior of the 
intensity spectrum, xdN^ /dx oc l/fc"*". Therefore, in this 
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FIG. 7: (Color online) Radiation intensity xdN^/dx for 
Ejct = 20 GeV, 100 GeV and 2000 GeV quark jets. Note 
tlie difference to first order in opacity in the large x = 
/p^ behavior for Bertsch-Gunion, initial-state and final- 
state bremsstrahlung. 



case, energy loss depends logarithmically on the jet en- 
ergy, as first derived in [ll, [l^ ■ 

The intensity spectrum, shown in Fig[7l is also impor- 
tant for understanding the effect of a heavy quark mass, 
Mq, on the magnitude of the energy loss, Ai^rad- Recall 
that the modulation of the 2D propagators and inter- 
ference phases is oc x'^M^, see Eq. ([7i|) and Ref. [1^. 
For the bremsstrahlung intensity, which has a negligible 
contribution in the region of large x where the mass cor- 
rection is significant, we expect that there will be no dif- 
ference between AE^ad for light and heavy quark. In con- 
trast, if the kinematically allowed region in a; = k'^ /p'^ 
contributes equally to the radiation intensity, as in the 
Bertsch-Gunion case, or the cancellation of the large x 
radiation is partial, as in the case of initial-state energy 
loss, then the dependence of Ai^rad on the heavy quark 
mass should remain. Numerical results for the three cases 
discussed in this paper are shown in Fig. [S] We used 
the same choice of parameters that describe cold nu- 
clear matter in Fig. ([7]). Quark masses for light, charm, 
and bottom quarks have been set to: Mu,d = GeV, 
Afc = 1.3 GeV, and Ah = 4.5 GeV. 

It is seen that, for final-state energy loss in cold nuclear 
matter, the charm quark AiSj-ad becomes comparable to 
the hght quark Ai;rad at E-^^t ~ (p'^/'^) = 20 GeV. For 
the much heavier bottom quarks, equality is reached at 
around 100 GeV. These results, given the crude energy 
binning used here to cover a very large dynamic range, 
5 GeV < p'^ /2 < 10® GeV, are not inconsistent with sim- 
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FIG. 9: (Color online) Differential in the gluon trans- 
verse momentum |k| bremsstrahlung multiplicity distribution 
dN^ /dyd^ln for lightcone momentum fractions x = /p^ = 
0.03,0.1 and 0.3. We show for comparison quark jets of 
Bjct = 100 GeV and 2000 GeV and the Bertsch-Gunion, 
initial-state and final-state spectra for comparison. 



FIG. 8: (Color online) Mass dependence of the fractional en- 
ergy loss AEiad/E versus the jet energy, Ejet- We studied 
massless quarks, Mq — GeV, charm quarks Afc = 1-3 GeV 
and bottom quarks Mb — 4.5 GeV. 



ilar fin ding s for heavy quarks traversing a quark-gluon 
plasma 23'|. On the other hand, for the Bertsch-Gunion 
case and the physically relevant case of initial-state en- 
ergy loss, the mass dependence persists at any jet energy, 
see Fig. HI 

The last numerical result, presented in this manuscript, 
is the fully differential gluon bremsstrahlung distribu- 
tion dN^ / dydP'k for select values of x — jp^ ver- 
sus the gluon transverse momentum. We have chosen 
X = 0.03, 0.1 and 0.3 to cover both small and large values 
of the gluon lightcone momentum fraction. Two differ- 
ent jet energies, Ejet = 100, 2000 GeV are also shown in 
Fig. [21 Massless quarks were used as an example in this 
study. The most characteristic feature of the Bertsch- 
Gunion energy loss is that at large values of the momen- 
tum, |k|, the bremsstrahlung spectrum behaves ^ 1/k'*. 
Note that this behavior is different when compared to 
the ^ 1/k^ for hard bremsstrahlung. The turnover in 
the growth of the spectrum at small |k| is controlled by 
the largest of the non-perturbative scales in the problem, 
in our case, Qq. Note that there is no real difference as 
a function of x or -Ejct, as expected from Eq. (jSip to first 
order in opacity. 

At the other extreme is the final-state medium-induced 
differential gluon distribution, see the right bottom panel 
of Fig. O The cancellation of the small |k| radiation of 
the basic Bertsch-Gunion term clearly becomes more ef- 
fective at larger values of x. Given the spectrum steeply 
falling in |k|, this cancellation is reflected in the large 
X suppression seen in Fig. [7] and the qualitatively dif- 



ferent behavior of the final-state Ai^rad when compared 
to QED or the Bertsch-Gunion limit. Initial-state en- 
ergy loss exhibits a cancellation which, however, is lim- 
ited at large values of x, see the middle bottom panel 
of Fig. [9l This explains the finite reduction of the en- 
ergy loss relative to the incoherent limit of the medium- 
induced bremsstrahlung in QCD. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper, in the framework of the reaction oper- 
ator approach 18'|, we derived the full solution for the 
medium-induced gluon bremsstrahlung for asymptotic, 
t — ztoo, on-shell fast partons (Bertsch-Gunion) as well 
as partons that undergo hard scattering [Q^ ^ Qq > 
/i^ > Aq^^ ) to produce high-pT or high-Er particles 
and jets (physically relevant case). The double differen- 
tial intensity spectrum, xdN^ / dxd^\i, was derived as an 
infinite opacity series: an expansion in the correlation be- 
tween the sequential multiple scatterings in the medium. 
Our general results are suitable for further analytic ap- 
proximations and/or numerical simulations. 

One of the main findings of this work is that, in con- 
trast to final-state energy loss, the cancellation of the 
Bertsch-Gunion radiation for initial-state Ai^rad due to 
the non-Abelian LPM effect is less effective. While there 
can be a large (in our examples six-fold) reduction of the 
total AE'rad relative to the incoherent Bertsch-Gunion 
case, such a cancellation does not qualitatively change 
the path length, c>c L, and the energy, (x E, dependence 
of the radiative energy loss for finite nuclei in the case 
of few interactions. It is clearly the incoherent Bertsch- 
Gunion regime that sets the upper limit for the amount 
of energy loss induced by strongly interacting matter. Of 
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the two cases relevant to hard probe physics, it is initial- 
state energy loss that by far dominates over final-state 
energy loss for the same medium parameters, such as the 
gluon mean free path, Ag, and the squared momentum 
transfer, /i^. 

The reason that final-state energy loss has been found 
to play a dominant role in A-f-A collisions is the large 
density of the QGP co-moving with the jet. Moreover, 
the effect of final-state AE'iad is amplified by the steeply 
falling spectra of outgoing partons, ^ ^/Px with n > 4 fli|. 
This observation implies that, in a region where the flux 
of incoming partons is rapidly changing as a function of 
the pt/Et of the final-state particle/jet, the cold nuclear 
matter AE^^ad can lead to a significant and observable 
suppression of hadron production [l^, US] • Such regions 
of phase space in p-l-A and A-l-A reactions are the ones 
near the projectile rapidity. In such kinematic domains in 
the rest frame of the target nucleus the incoming partons 
are of almost asymptotically high energy. The knowl- 
edgeable reader can easily verify that for pT — l — 2 GeV 
at rapidity y = 4, a region accessible by the STAR col- 
laboration at RHIC 7], the incoming partons are of en- 
ergy E > 20 TeV. Therefore, the contribution to AE 
from final-state radiative energy loss and coUisional en- 
ergy loss is completely negligible. In fact, due to the ra- 
pidity boost, the incoming and outgoing parton energies 
are very large everywhere, except in the backward region 
near the target rapidity. Theoretical results, derived in 
this manuscript, can be used to develop a more complete 
and consistent pQCD phenomenology of proton-nucleus 
collisions [l3|- Given the short time scale of dynamical 
cold nuclear matter effects, tco\ «C icq ^ L^^^ , these 
must also be incorporated in the description of hard pro- 
cesses in A-l-A reactions. 

We also argue that a new application of cold nuclear 
matter-induced radiative energy loss can be related to 
particle production at backward rapidities. It has been 
established, through the measurements of the PHENIX 
collaboration 25] at RHIC, that the forward rapidity 
hadron suppression in p-|-A reactions is correlated with 
a backward rapidity enhancement. There are no mod- 
els at present that can consistently account for this ef- 
fect. One possibility is that such enhancement in the in- 
clusive particle yields may arise from the induced gluon 
bremsstrahlung. What makes such a scenario plausible is 
the fact that the gluon yield, dN^ /dxcPk, is dominated 
by small x radiation and the destructive LPM effect can- 
cels the small |k| part of the induced gluon spectrum. 
For the final-state AE'rad, this contribution is expected to 
yield an enhancement of the soft, large- angle hadrons as- 
sociated with an away-side triggered jet [2a]. For initial- 
state cold nuclear matter energy loss, this corresponds 
to gluons emitted preferentially in the backward rapidity 
(or near target rapidity) region, see Eq. (|76p . Since at 
moderate /cy the power behavior of the gluon spectrum 
is similar to that of hard scattering, oc it may have 

a sizable contribution to the inclusive particle yields. 

We emphasize that the results derived in this pa- 



per are applicable to both cold nuclear matter and the 
quark-gluon plasma. However, the physics situation of 
asymptotic — cx), +oo jets propagating through the QGP 
(Bcrtsch-Gunion) or jets undergoing hard scattering at a 
finite time t — L, having penetrated the QGP, (initial- 
state) cannot be experimentally realized at present. Nev- 
ertheless, the formal radiative energy loss solutions ob- 
tained here are general enough to describe both cases of 
interest. 

In summary, we anticipate that the results derived in 
this paper, in particular the ones related to initial-state 
cold nuclear matter radiative energy loss, will play an 
important role in a consistent many-body pQCD descrip- 
tion of hard processes in high energy reactions of heavy 
nuclei and complement the well-developed theory and 
phenomenology of final-state energy loss in the quark- 
gluon plasma. 
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APPENDIX A: SIMPLIFYING THE ENERGY 
LOSS CALCULATION FOR ENERGETIC 
PARTONS 

In this paper we adopt the notation: 

p^ = b+,p-,p]= [pO+p-,/-p-,p] , (Al) 

for the lightcone momentum components. The scalar 
product between two 4-vectors is then given by: 

p-q^^{p^q^ +p^q^) -p-ti- (A2) 

The soft gluon approximation that we use implies: 

fc+ < p+ (A3) 

for the positive lightcone momentum of the gluon rela- 
tive to the lightcone momentum of the jet. Evidently, 
the gluon multiplicity should be dominated by small 
X — /p'^ gluons, which we can verify a posteriori in all 
cases. The gluon intensity spectrum, however, can have 
contributions from the large x region of phase space. For 
the transverse momenta we require: 

^^-^^"^^^ « k+ . (A4) 

This constraint implies that at the emission vertex the 
transverse momentum k — of the gluon is small 
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and the large transverse momenta k are accumulated via 
multiple scattering, which can be large angle. Finally, 



p+ 



0< 



(A5) 



and such terms are neglected. The polarization vector 
for the physical final-state gluon is given by 



e{k) 



(A6) 



so that e{k)-k — 0. With n momentum transfers from the 
medium and the approximations outlined in Eqs. (jA3|) - 
(jA5p . the kinematic part of the gluon emission vertex 
reads: 



2p+(fc+)"-^ex-(k-q., 



(A7) 



The color factor, associated with Eq. (|A7|) . is 

= di is the color matrix at the emission vertex. Thus, 
a common factor —2igse±^ ■(■■■) can be factored out from 
all amplitudes. 



While the integrals over the longitudinal position of the 
scattering centers, Zi, have to be taken explicitly, a major 
simplification at the level of squared amplitudes occurs 
when we consider the average over the corresponding im- 
pact parameters = — Xq. This average is done in the 
local rest frame of the medium containing colored scat- 
tering centers. For a scattering center, i, that appears in 
a direct interaction the impact parameter average takes 
the form: 



A, 



A;(0, qi)e 



(2^) 



(27r)2 (27r)2 Ai_ 



1 dop 



(J el d^qi 



(A8) 



which constrains the transverse momentum exchanges q.; = q^ to be equal in the amplitude and its conjugate. For a 
double-Born interaction: 



d^h,. 



Ai 



(-0 



■(-1) 



(27r)2 



Ai(0,q^)e-*'i-'^'(-0 



d^ 



Ai(0,q,)e 



(27r)2 (27r)2 



A, 



Ael{0,q^)Ael{0,q'^) 



d'^q^ 

(2^)2 



1 da. 



(J el rf^q* 



j d^q^6''{ci, + c{,) 



(A9) 



The momentum transfers are constrained here to be equal and opposite. In Eq. (IA9P we have used ^e;(0,qi) 
^g,(0,qi) and Ai(0,qi) ^ AeiiO, -qi)- ^ 



APPENDIX B: AMPLITUDE ITERATION 
TECHNIQUE TO SECOND ORDER IN OPACITY 
FOR ON-SHELL PARTONS 

We illustrate the iteration of gluon emission am- 
plitudes to second order in opacity and calculate the 
Bertsch-Gunion case for direct comparison to the gen- 
eral gluon emission result for soft scatterings Eq. ([5T|) . 
One should first recall that for the Bertsch-Gunion case 

Go - , (Bl) 

since the asymptotic on-shell quarks do not radiate with- 
out the medium-induced acceleration. For the two rank 
1 classes, we apply once the direct and virtual insertion 
operators, Di from Eq. pT|) and Vi from Eq. (P?)) . to the 



incoming parton, Eq. (|Bip . and obtain: 

Z?iGo = -Bie'"°-Mc,ai] , (B2) 

14 Go = -^Bie^-«-ic . (B3) 

We note that the results are simple since only the new 
color current Bi contributes in the absence of a non-zero 
initial gluon amplitude. 

To second order in opacity, we build upon the am- 
plitudes from Eqs. (|B2[) and (jB3|) . Some of the rank 2 
classes are obtained from rank 1 through relabeling, i.e. 
i^2Go = i^iGo(l ^ 2), V2G0 = V^iGo(l ^ 2). The rest 
are readily derived from Eqs. (jB2p and (jB3p through our 
iteration scheme Eqs. P7|) and (|27p : 
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V2D1G0 



-B2(i2) e 

Cr + Ca 



i(<^o.2-<.2.2i)[[c^^2],ai] , (B4) 
B,e'-^°'' [c,ai] 



Ca 



To second order in opacity we have: 

C{ = a2aic , 63 = a2 [c, ai\ 
Cg = [c, fli , 64= ca2ai 



B2 e'"«^=^ cai 



,(2) 



[c, 02] , aij 



(B13) 



+B2(i2) e'^'^o'^^ ""^"^^i) a2[[c, 02], ai] , (B5) ^he traces over the color factors yield: 



D2V1G0 = -^Bie^'^''^^a2C+^B2e'"«"Mc,a2] 



'^^2(12) e 



i(iA;oZ2-1^2Z2l) 



C, a2j 



(B6) 



r2 



B2 e'"°^== c 



-r ^2(12)6 



i(aJ0Z2-I^2Z2l) 



(B7) 



With this explicit construction of the relevant classes of 
diagrams, we can compute the differential gluon emission 
up to second order in the opacity expansion. 

A necessary side step in the brute force approach is 
the evaluation of the color factors using the techniques 
in Refs. [13, HI]. We denote by Cr the quadratic Casimir 
of the representation of the incident parton. For SU(iVc), 
following the standard normalization for the generators 
we have 

- 1 

Cf = Ca = N,. (B8) 

We recall that in our notation Tc = c for brevity. In the 
absence of interactions we have only one color matrix 
associated with the gluon emission vertex from a hard 
scatter, 



with 



TtC["^^c["^ =CrDj 



(B9) 



(BIO) 

To first order in opacity we have to consider the following 
color factors 

C[^'^^cai, 4^) = [c,ai], = aiC , (Bll) 



with 



Tr C['^k['^ - Cf,DR , C^^^'C^^^ = CrCaDr , 



(i)t^(i) 



T,c[')tci'^ = ^CRDR, 



TtC['^^4'^^ [Cr--Ca]CrDr, 



1 



(B12) 



TrCf^M^) =Tr cf^cf^C^i^^, 
TrC^^^^C^^^ = Tr Cf^^cf^ = CaC^Dr , 



TrC^^^^C^^^ = C\CrDr , 



Tr cf^^C^^' = -Tr C^'^'C 



(2)t/.(2) 



3 "-4 



-CACj^DR , 



TrC}''^4'^ = -TrC^^^Tc(2 



= --^ yCR - -C'a ] C'aGrDr 

TrCf ^^Cf ^ = (cr - ^Ca^ CrDr , 
Tr Cf^^cf^ = , 



1 



-CaCrDj 



(B14) 



Although some of the color factors, such as Cg^'' and 
do not appear in the amplitudes listed above, these will 
prove useful in the calculation of the realistic case of in- 
coming partons that undergo hard scattering to produce 
final-state jets. 

To first order in opacity, the radiation from either 
quark or gluon jets reads: 



^+nI'\bg) ^ 1 



oc 



dfc+d^k Dr 



{ Tr [iD,Go)^D,Go 
+ {{Go)^ViGo + h.c.)] ) .(B15) 



In Eq. (|B15p . (■ • ■) denotes the average over the momen- 
tum transfers from the medium. Including the phase 
space factor, the strong coupling constant ag and inte- 
grating over the longitudinal position of the scattering 
centers we find: 

, + ivW(BG) _ CRa, /-^ dAzi f 2 



dk+d^k 



Jo >^g{zi) 



d'qi 



x^—^[\B,\'] . (B16) 



(Tei(zi) d^qi 

Here, we recover the Bertsch-Gunion result to one scat- 
tering center [Til] and show that it is equivalent to first 
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order in opacity. One important difference between the 
initial- and/or final-state radiation and on-shell jets is 
that in the latter cases the first order in opacity does not 
lead to any interference/coherence effects. 

The first non-trivial correction to the Bertsch-Gunion 
result can also be directly calculated. Using the ampli- 
tudes, Eqs. (|B2|) - (|B7|) . and the calculated color factors, 
Eqs. (|BTT|) - (|BT4l) . we find: 



1 



+ {{D2GoyD2ViGo + h.c.) 

+ {{DiGoyV2DiGo + h.c.) 

+ {{Go)^V2ViGo+h.c.) 

+ mGo)^ViGo+h.c.)] ) .(B17) 

The final result for the first non-trivial coherence correc- 
tion reads: 



n^'\bg) 

dk+d^k 



C 



dAzi 



1 daei 



L-Azi 



dAz2 
^3(^2) 

1 dffei 



aei{zi) d'^qi J (Jei{z2)d^q2 

[ |B2(i2) P - |Bi |2 + 2B2 ■ B2(12) C0S(W2221) ] • 

(B18) 



In a medium with no sharp boundaries we can extend the 
integration limits as follows, z(0) — > —00, z{L) +00. 
The disappearance of the interactions is guaranteed by 
\g{zi) ^ 00 as p{zi) — > 0. Direct comparison to the 
full recurrence solution, Eq. (|5ip . can now be made with 
consistent results, as expected. 



APPENDIX C: AMPLITUDE ITERATION 
TECHNIQUE TO SECOND ORDER IN OPACITY 
FOR INITIAL-STATE RADIATION 

The case of initial-state energy loss, as emphasized 
in the general derivation section, differs from the full 
Bcrtsch-Gunion solution by the presence of the additional 
hard scatter at position L. While the t = —00 state for 
the incoming on-shell parton remains the same, the hard 
acceleration yields an additional term: 

H = ( - 2 J He^""'" cTel{A^,...^J . (CI) 



In the absence of soft interactions. 



Go + H = He' 



(C2) 



The dependence of H on the sequence of interactions 
preceding the hard scatter is implicit on our notation. 
If we want to explicitly verify the general result to sec- 
ond order in opacity, we need to calculate only a few 
additional terms relative to Appendix [B) The relevant 
amplitudes are modified as follows, 

DiGo + H = -Bi e"^"'' [c, ai] 4^ Ue"^°'^ cai , (C3) 



Ca 



ViGa + H = — Bi e*"°^i c - ^He 



(C4) 



and 



D2D1G0 + H = -Bie^'^''^ia2[c,ai] 
-Bae*"""^ [c,a2]ai 
-B2(i2) e*('^''^^""^^^i) 

Cr + Ga 



[c, a2\,ai\ 



(C5) 



V2D1G0 + H = 



Bie^'^o^i [c,ai] 



Ga 



B2 e"^°'--' cai 



+B2(i2) e 



2 

D2V1G0 + H = -^Bie^"°^ia2C 



02 c, a2\,ai\ 



(C6) 



R 



G 

2 

Ga 



B2 e'^"'' [0,02] 



B2(i2) e^("«^^-"^^-)[c,a2] 



V2V1G0 + H 



-fHe--ca2 
Ga{Gr, + Ga) 



(C7) 



Bi e''^°^i c 



GrGa 



B2 e" 



^2 

4 



-02(12) e 



(-12 

4 



(C8) 



in comparison to Eqs. (jB2[) - (jB7[) . Having also calculated 
the color factors for this case in Appendix [Bl we obtain 
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fc+^|+^ cx —{TT[{D2DiGo + H)^D2DiG„ + H) + {{D2Go + HyiD2ViGo + H)+h.c.) 

+ {{DiGo + H)\V2DiGo + H)+ h.c.) + ((Go + Hy{V2ViGo + H) + h.c.) 
+ {{V2Go + H)^{ViGo + H)+h.c.)]) . (C9) 

I 

Evaluating the additional terms, arising from the inter- pie Bertsch-Gunion sources, and showing explicitly that 
ference of the hard gluon bremsstrahlung with the multi- the term oc cancels, we obtain: 



^"^^ - ^ / ^^^^^^ [IBiP -2H.B,cos(.o..r)] • (CIO) 



>2(12) 



dk+d'^k 7r2 \g{zi) Jq Ag(z2) 7 aei{zi) d'^cii J aei{z2) d'^q2 

-|Bip + 2B2 • B2(12) COs(a'2Z2l) - 2H • B2(12) COS(WOZL2 + i^2Z2l) + 2H • Bi cos{uqZli) ] ■ (Cll) 



I 

These results coincide with the general opacity expan- bremsstrahlung, see Eq. (|63p . 
sion series for initial-state medium-induced non-Abelian 
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